Evolution of Fluctuation in relativistic heavy-ion collisions 
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We have studied the time evolution of the fluctuations in the net baryon number for different 
initial conditions and space time evolution scenarios. We observe that the fluctuations at the freeze- 
out depend crucially on the equation of state (EOS) of the system and for realistic EOS the initial 
fluctuation is substantially dissipated at the freeze-out stage. At SPS energies the fluctuations in 
net baryon number at the freeze-out stage for quark gluon plasma and hadronic initial state is close 
to the Poissonian noise for ideal as well as for EOS obtained by including heavier hadronic degrees of 
freedom. For EOS obtained from the parametrization of lattice QCD results the fluctuation is larger 
than Poissonian noise. It is also observed that at RHIC energies the fluctuations at the freeze-out 
point deviates from the Poissonian noise for ideal as well as realistic equation of state, indicating 
presence of dynamical fluctuations. 

PACS numbers: 25.75.-q,05.40.-a,12.38.Mh 



I. INTRODUCTION 



Heavy-ion collisions at relativistic energies offer a 
unique environment for the creation and study of 
the Quark-Gluon Plasma (QGP) phase of nuclear 
matter. In the colhsion process the system may 
go through a phase transition whereby it evolves 
from a very hot and dense QGP state to normal 
hadronic matter. A characteristic feature of this 
process is that the system experiences large event- 
by-event fluctuations in thermodynamic quantities 
such as temperature, number density etc. Our abil- 
ity to observe a characteristic fluctuation in various 
observables, in present day experiments, has been fa- 
cilitated by the large number of particles produced in 
the relativistic heavy-ion collisions at the Super Pro- 
ton Synchrotron (SPS) and the Relativistic Heavy- 
Ion Collider (RHIC) along with the advent of large 
acceptance detectors at these experiments [Q. The 
near Gaussian distributions of the experimental ob- 
servables like multiplicity, transverse energy 
transverse momentum and particle ratios H have 
provided an opportunity to relate these quantities to 
thermodynamical properties of matter, such as spe- 
cific heat 0] and matter compressibility Fluc- 
tuations as a signature of quark-hadron phase tran- 
sition has been experimentally studied extensively 
but no definite signature has been observed in this 
sector so far |, |, 

The main task for using fiuctuation as a probe 
of QGP phase transition would be to identify an ob- 
servable whose fiuctuation will survive from the time 
of formation of plasma till they get detected after the 
freeze-out. In this context, it was suggested that the 
study of the fluctuation in conserved quantities e.g. 
net baryon number, net charge or net strangeness 
will be very useful |, |, The fluctuation in 



any conserved quantity (O) is given by p 
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where, /i is the associated chemical potential and 
T is the temperature of the system under consider- 
ation. Here (O) represents the mean value of the 
conserved quantity over a large number of events 
and (AO) denotes the deviation of its value from the 
mean, event-by-event. It has been shown that the 
fluctuation in these conserved quantities with and 
without QGP formation are very different. Hence, 
study of fluctuations in these conserved quantities 
can be a good signature of transition from confined 
to de-confined state of matter. 

However, in the earlier calculations |l^ ideal 
gas EOS for both the QGP and hadronic gas was 
considered. For the QGP initial state, the evolu- 
tion of fiuctuation during mixed phase have also not 
been taken into account. Finally, modifications of 
hadronic spectral functions in the thermal bath was 
ignored. 

The aim of this paper is to study the evolution of 
the fluctuation in conserved quantities from time of 
formation to the time of detection. In addition to the 
QGP and hadronic initial state we consider an im- 
portant case where the hadrons are formed at the ini- 
tial stage but the spectral function of the hadrons are 
different from their vacuum properties. We also dis- 
cuss the sensitivity of the results on different EOS. 
Here we concentrate only on the net baryon number 
fluctuation and the results are presented in terms of 
the experimentally measured quantity, {ANb)"^ /Nt, 
where Ni, is the number of baryons and ANf, is the 
fluctuation in the net baryon number. 

The paper is organized as follows. In the next 
section we discuss the possible evolution scenarios. 
Then the initial conditions used for solving the evo- 
lution equations for various evolution scenarios are 
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presented in section III. In section IV we study the 
evolution of the system. Section V and VI is devoted 
for the results at the freeze-out for SPS and RHIC 
energies respectively. And finally in section VII we 
present the summary and conclusions. 



III. INITIAL CONDITIONS 

The initial conditions in terms of initial tempera- 
ture {Ti) can be set for the three scenarios from the 
following relation: 



II. EVOLUTION SCENARIOS 

We consider the following three possible evolution 
scenarios. 

1. QGP scenario: The system starts evolving 
from an initial QGP state formed at a time 
Ti and temperature Ti. It then expands, hence 
cools and reaches the critical temperature Tc 
at a time Tq. In the mixed phase the cool- 
ing due to expansion is compensated by the 
heating of the system due to the liberation of 
the latent heat in a first order phase transi- 
tion. Hence the temperature remains constant 
at Tc (super-cooling is neglected here). Af- 
ter the mixed phase ends at a time r/j, further 
expansion takes place and finally the system 
disassembles at a time Tf and temperature Tf. 
At this stage, called the freeze-out stage, the 
mean free path of the particles are too large 
to have any further interactions, and those are 
detected experimentally. 

2. Hadron gas scenario: The hot and dense sys- 
tem is formed in the hadronic state at a time 
Ti and temperature Ti and the system expands 
from the initial to freeze-out state (at time Tf 
when the temperature is Tf) without a phase 
transition. 

3. Hadron gas with mass variation: The other 
possibility is that the system may form in the 
hadronic state as in (2) above, but the spec- 
tral function of the hadrons are different from 
their vacuum counterparts. Among others, in 
the present case we will consider the shift of 
the pole of the hadronic spectral function ac- 
cording to the universal scaling law proposed 
by Brown and Rho |l3| : 



dS = 



mil 



T2 



(2) 



where m*^ (rrih) is the in-medium (vacuum) 
mass of the hadrons (except pseudo-scalar) 
and the index A takes a value between and 
1. Here we choose A = 1/6 according to the 
well known Brown- Rho scaling fl^. 



(3) 



where dS{dN) is the entropy (number) contained 
within a volume element AV — ir Ti di], i? as the 
radius of the colliding nuclei, g^f / is the effective sta- 
tistical degeneracy, C(3) denotes the Riemann zeta 
function and 77 is the space-time rapidity. For mass- 
less bosons (fermions) the ratio of dS to dN is given 
by 27r4/(45C(3)) 3.6 (4.2), which is a crude ap- 
proximation for heavy particles. For example the 
above ratio is 3.6 (7.5) for 140 MeV pions (938 MeV 
protons) at a temperature of 200 MeV. 

First we consider the situation at SPS energies 
where dN/dy ~ 700 for Pb+Pb collisions. For above 
three scenarios taking g^f / as given in Table ^ we 
obtain the initial temperature by using Eqn. ^j. The 
values of the initial temperatures are given in Ta- 
ble | 

The chemical potential at the initial state is fixed 
by constraining the specific entropy (entropy per 
baryon) to the value obtained from the analysis of 
experimental data. The specific entropy at SPS is 
about 40 for Pb + Pb collisions |l|, |l|, 0. 
The net baryon number can be calculated using the 
baryon (anti-baryon) number density, rib (nf,) given 

by, 



{2nf 



(4) 



where g is the baryonic degeneracy. We take g = 4 
for proton and neutron. f(p) is the well known 
Fermi-Dirac distribution. 



fip}= [exp(^{E±^i)/T)+^ 



(5) 



where ^{—^) is the chemical potential for baryon 
(anti-baryon), E = -^/p^ + m^ and = rn,{iJL ^ 

For the QGP scenario we take the mass of the 
quarks as, niq^ = m^c^ + 'mgth'^, where uigth is the 
thermal mass and mqc is the current quark mass. 
We have taken vacuum mass and effective mass for 
the hadrons (given by Eqn. |2|) for cases 2 and 3 re- 
spectively. We obtain the value of the chemical po- 
tential to be 132 MeV, 340 MeV and 105 MeV for 
cases (1), (2) and (3) respectively. These are also 
summarized in Table 0. 

Having fixed the initial temperature and chemical 
potential we now calculate the initial fluctuations 
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TABLE I: Initial conditions and the initial values of fluc- 
tuation for the three scenarios. 



Initial Values/Scenarios 


1 


2 


3 


9eff 


37 


15 


24 


T, (MeV) 


196 


264 


226 


^^^ (MeV) 


132 


340 


105 


{ANtiT,))ys 


0.014 


0.029 


0.061 


{ANtiT,))''/Nt,y 


0.56 


1.16 


2.46 



for the three scenarios. 

1. Quark Gluon Plasma: From the initial net baryon 
number, one can easily calculate the net baryon 
number fluctuations using Eqn. The fluctuation 
in the net baryon number in the QGP phase at time 
Ti, temperature Ti and chemical potential ^i, can be 
shown to be [noi: 



(A7V,(r,))2 



QGP 



2V o 

rpo 

9 ' 



1 



]_fJM_ 
3 VttT, 



(6) 



With the initial conditions as discussed above, the 
initial fluctuations for the QGP scenario turns out 
to be (A7V;,(Ti))Qqp = 35. It may be noted that the 
value of g in Eqn.]j is taken to be 12, for two flavor 
case. 

The entropy density in the QGP phase is calcu- 
lated from Eqn. |^ with g^^ff — 37. The total entropy 
of the system is given by: 



S = V- 



90 



2500. 



(7) 



For fixed initial entropy of the system the initial 
temperature is different for the three cases because 
of the different values of the geff- In the present 
work we solve the evolution equation for ideal 
fluid neglecting entropy generations due to various 
viscous effects. The total entropy is kept constant 
for all the three scenarios as it is obtained from 
the number of particles per unit rapidity measured 
experimentally. So the fluctuation in net baryon 
number per entropy is {ANb{Ti))QQp/ S — 0.014, a 
value similar to that obtained in [ |lO| . In this case 
the fluctuation per unit baryon, {ANi,{Ti))QQp / N^^y 
= 0.56 where N^^y is the net baryon number per 
unit rapidity, dNb/dy ~ 62 for SPS energies. 

2. Hadron Gas : The net baryon number fluctuation 
in the hadronic gas can be calculated using Eqn. I 
and is given by: 



gV \2mTi 

T 



3/2 



exp{—m/Ti) cosh{^/Ti] 



(8) 



Substituting the values of Ti, /ii and m = 938 MeV 
for the nucleons, we obtain the fluctuation in net 
baryon number, (AA^f,(ri))^Q ~ 72. The ratio of 
the fluctuation to total entropy is {ANb{Ti))^Q / S = 
0.029. This value is also similar to the one obtained 
in [|l0). In this case the value of {ANb{n))'^Q/ Nb^y 
= 1.16. 

3. Hadron gas with mass variation in medium: The 
initial fluctuation for this case can be obtained using 
the Eqn's ^ and H together with Eqn. || used for the 
value of mass. The initial conditions of Ti — 226 
MeV and fii — 105 MeV are constrained to repro- 
duce the hadronic multiplicity and the entropy per 
baryons. In this case we get (AA^;,(ri))^Q j^^. =153 
and {ANbiTi))lQ^,/S = 0.061. The value of 
(A7Vfc(T,))|G,^./iVb,, =2.46 here. 

We observe that for the baryon number fluctua- 
tions corresponding to the above three scenarios: 



(ANbin)) 
{ANbin))^ 



2 

HG 



2 and 



(ANbin))^ 



HG,m* 



QGP 



{ANbir,))' 



4. (9) 



QGP 



This clearly shows that in the initial stage, there is a 
clear distinction between the three cases. The initial 
values of fluctuations are summarized in Table |. It 
will be of interest now to see, if these fluctuations 
(and the differences given by Eqn. |^) survive till the 
freeze-out. Next we discuss the evolution of this ini- 
tial fluctuation for the three scenarios. 



IV. EVOLUTION OF THE INITIAL 
FLUCTUATION 

We follow ref. [|l^ to study the proper time (r) 
evolution of the baryon fluctuation in the space-time 
rapidity interval A rj for different EOS as well as ini- 
tial states as mentioned above. At the freeze-out 
point the fluctuation measured experimentally con- 
tain the residue of the initial fluctuation which sur- 
vived the space time evolution and the fluctuations 
generated due to the exchange of baryons with the 
adjacent sub- volumes. We discuss each of the above 
cases separately in the following sub-sections. 



A. Dissipation of the initial fluctuation 

The difference in baryon flux (A n v) originating 
from different densities inside and outside of the sub- 
volume {At A ri) leads to the following differential 
equation for AA^f,(= AnvAr): 



dANb 
dr 



V ANb 
2A^ r 



(10) 
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FIG. 2: The variation of energy density as a function of 
temperature. The result shown by solid line is obtained 
from Eqn. Open circles show the lattice results with 
out the error bars. 



FIG. 1; Variation of average thermal velocity with tem- 
perature for three different scenarios. The results are 
fitted by nth-order polynomials. (1) For QGP scenario 
the fitting parameters are ao = 0.93, a\ — 0.11, 
a2 — — 0.56, as = 1.32 and 04 — 1.20. 
(2) For hadronic gas scenario with fit parameters as 
ao = 0.19, ai = 4.0, 02 = - 14.4, 03 = 32.2 
and a4, — — 31.2. (3) For hadronic gas with mass vari- 
ation in medium with fit parameters as ai = — 4889.8, 
02 = 76879, as ^ - 0.640E + 06, 04 = 0.297S + 07, 
as = - 0.73.E + 07 and ae = 0.74E + 07. 



where v is the average thermal velocity of the parti- 
cles under consideration. 

The solution to the differential equation |l^ from 
initial time t; to final time tj is then given as: 



A7Vb(r/) = AiVb(Ti)exp 



1 



dr 



2A7J T 



w(r) 



(11) 

The average thermal velocity at a given tempera- 
ture can be calculated by using the following equa- 
tion, 



Jd^pfip) 



(12) 



The values of w as a function of temperature for the 
three cases are shown in Fig. |l|. The curves are 
fitted by the polynomials of the form ap -I- aiT -f 
aaT^ + asT^ + a^T^ + a^T^ + a&T^. The values of 
the parameters are shown in the caption of Fig. |^. 
The average velocity in QGP is found to be sub- 
stantially larger than the velocity in hadronic gas. 



However, in case of the hadronic system where the 
effective masses of the baryons (neutron and proton 
here) approach zero at Tc, their velocities approach 
the velocity of light. 

The equation of state plays a vital role in decid- 
ing how the fluctuations in the net baryon number 
will evolve. The variation of energy density (e) as a 
function of temperature is obtained from the 3-flavor 
lattice QCD results which is parametrized as fol- 
lows: 



T^A tanh{B{^f), 

^ C 



(13) 



where the values of the parameters. A, B and C are 
12.44, 0.517 and 10.04 respectively. Note that the 
effect of net baryons in the EOS is neglected here. 
Fig. ^ shows the variation of the energy density with 
temperature. The increase in the effective degener- 
acy near Tc can be obtained from the hadronic phase 
with effective masses varying with temperature as in 
Eqn. I The increase in the effective degeneracy orig- 
inates from the heavier hadrons going to a massless 
situation ( see also ref. ^). 

The evolution of the system under the boost in- 
variance along the longitudinal direction is governed 
by the equation 123], 



rf7 



P 



0; P 



(14) 



where Cs is the velocity of the sound in the medium. 
To perform the integration in Eqn. |ll|we need dr/r, 
which can be obtained as 



dT 
T 



adx 
xdT 



dT ■ 



-f{T)dT 



(15) 
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where, a = yq^, /3 = 4a and 



X 



tanhBiTjTcf 



(16) 



The expression for the second term within the square 
bracket of Eqn.nS is given by 



XdT 



1 



(17) 
//i-i 



Using the relation = y ^ = [3 + 
we get the velocity of sound corresponding to the 
lattice QCD calculations, 



1 



tJ cosh{B{^)C) sinh{B{f)C) 

(18) 

Now the fluctuation at the freeze-out point can be 
written as, 



AiVb(T/) = AiVb(r,)exp 



1 



2At] 



Tl 



fiT)dTv{T)j 
(19) 

In order to study the effect of the equation of state, 
we will present results for three different values of 



Cs'- (a) — 1/3 corresponds to the ideal gas case, 
(b) = 0.18 corresponding to an EOS of hadronic 
gas where particles of mass upto 2.5 GeV has been 
taken into account from particle data book and (c) 
cz obtained from Eqn. |l|. It may be mentioned that 
for the ideal gas the rate of cooling is faster than the 
other two cases, (b) and (c). 



Now we shall calculate the total dissipation of 
fluctuation at the freeze-out temperature. We 
shall consider the freeze-out temperature to be 120 
MeV |lj, HJ] and a critical temperature for the 
QGP transition to be 170 MeV In all the 

three cases below we have taken the value of Arj = 1. 



1. Quark Gluon Plasma: For the QGP initial state, 
the dissipation equation has three parts. In the first 
part, we calculate the dissipation from the initial 
temperature Ti to Tc (for the QGP phase). The sec- 
ond part is the dissipation during the mixed phase 
and in the final phase it is the dissipation from Tc 
(end of mixed phase) to the freeze-out temperature 
Tf. Hence the complete dissipation equation be- 
comes. 



A7Vb(r/) = AArb(r,)exp 



U ' v,gp{T)f{T)dT + VmUTc) ln{r) ^ jj vn{T)f{T)dT\\ (20) 



r 



where r is the ratio {gef f)qgp/gef f)h ~ 2.5 and 

The fluctuation in the net baryon number is 
evaluated for: (a) — 1/3, for which we 
get {ANb{Tf))%^p ^ 6, (b) cl = 0.18, giving 
{ANb{Tf))l^Q^ - 2.23 and (c) cl as given by Eqn.|l|, 
resulting in (AiVf,(T/))QQp ~ 0.54. 

The initial and flnal fluctuation values along with 
the evolution process of the values for the QGP 
scenario is shown in Fig. ^. Only the results for 

= 1/3 and as given by Eqn. ^ are shown in 
the figures for the clarity of presentation. We have 
checked that the results for = 0.18 lies between 
the above two cases. 



2. Hadron Gas 
have. 



In case of hadronic initial state we 



(a) For = 1/3, the contribution from the expo- 
nential is ~ 0.5. Hence the fluctuation in net baryon 
number is (A7V6(t/))|q is 18, (b) for = 0.18, we 
get (AiVf,(r/))yQ is 14.6 and if we take (c) from 
Eqn. ^ the contribution from the exponential term 
is ^ 0.06, a factor of 8 smaller than the value for 
case (a) . Consequently the fluctuation in net baryon 
number is {ANb{Tf))^Q is 0.22. The evolution of 
fluctuation for this case is depicted in Fig. ||. 



3. Hadron gas with mass variation in medium: For 
the hadronic initial state with mass variation. 



ANb{Tf) = A7Vfc(T,)exp 



^X^'<(T)/(T)dTj22) 



(a) For — 1/3, the contribution from the ex- 
ponential is ~ 0.53. Hence the fluctuation in 
net baryon number is (AAfb(T/))^. 



Am{Tf) = ANi,{T,)e^p 



2Ar] 



- 44, (b) for 

\ ^-^^^ J^^biTf Jf^, - 37 and (c) 

Vh{T)f(T)dT\2l f^^^''^S ''s from Eqn. |l8|, the contribution from the 
'j, / exponential term is ~ 0.08 a factor of 7 lower than 
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FIG. 3: The dissipation of net baryon number fluctua- 
tion for the different scenarios as a function of temper- 
ature. Note that temperature is plotted in decreasing 
order to reflect the evolution in time. (1) Variation of 
(AiVi,)^ for QGP scenario, (2) Variation of {ANtf for 
hadronic eas scenario and (3) Variation of (ANb)'^ for 
hadronic gas with mass variation scenario. The dashed 
lines correspond to results obtained for = 1/3, while 
the solid lines show the results obtained using the value 
of given in Eqn. |l^. 



the case (a). So the fluctuation in net baryon 
number is (AA'h(T/))^Q j^^, is 0.92. The evolution 
of fluctuation for this case is shown in Fig. 0. 



B. Generation of fluctuation with time 



{AtAt]) between time Ti and r/ is given by Eoj 



NriTf) = v{T)f{T)dTj (23) 



1. Quark Gluon Plasma: As in the previous case 
the fluctuations has to be evaluated in three parts 
in a QGP formation scenario. In the first part, we 
calculate the generation from the initial temperature 
Ti to Tc (for the QGP phase) . The second part is the 



150 




300 275 250 225 200 175 150 125 

Temperature (MeV) 



FIG. 4: The generation of net baryon number fluctua- 
tion for the different scenarios as a function of tempera- 
ture. Notations are same as that of Fig. |^. 



The baryon fluxes exchanged with neighboring 
sub- volumes can lead to generation of fluctuations, 
and is the main source that is to be detected exper- 
imentally at the freeze-out point. The total number 
of baryons leaving or entering {N^^) the sub- volume 



generation during the mixed phase and in the final 
phase it is the generation from Tc (end of mixed 
phase) to the freeze-out temperature Tf. Hence the 
complete evolution equation becomes, 



{AN,{Tf)f = Nt{T,) [^\^J^' Vqgp{T)f{T)dT + v^„{Tc) ln{r) + jj Vh{T)J{T)dTyj (24) 



r 



The fluctuation in the net baryon number is (AiVh(T/))QQp is 56.5, (b) = 0.18, resulting in 
evaluated for: (a) — 1/3, for which we get 
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(AAfb(T/))^QP is 60.4 and (c) from Eqn. |l| giv- 
ing iANt{Tf))lQp is 129. 

The initial and final values of the fluctuations 
along with its time evolution for the QGP scenario 
are shown in Fig. ^. 

2. Hadron Gas : In case of hadronic initial state we 
have, 



(a) For cl ~ 1/3, the fluctuation in net baryon 
number is (AiVb(T/))^(3 is 44. (b) For = 0.18, 
we get_(AA^b(T/))^Q is 50. (c) If we take from 
Eqn. |l^, the fluctuation in net baryon number is 
(AA^b(T/))yQ is 179. The evolution of fluctuation 
for this case is displayed in Fig. ^. 



3. Hadron gas with mass variation in medium: For 
the hadronic initial state with mass variation, 



TABLE II: {ANt f/Nb,y at freeze-out for three different 
evolution scenarios at SPS and RHIC energies. 



EOS /Scenarios 



ci = 1/3 To LO O 

SPS = 0.18 1.0 1.0 1.3 

Lattice 2.0 2.8 2.5 



= 1/3 



1.46 



RHIC cl = 0.18 1.96 
Lattice 3.17 



c2 = 1/3, iAN,iTf))l^/m^y ^ 1.0 (b) for - 
0.18, (ANb{Tf))^Q/Nby - 1.0 (c) Taking from 
^ - 2.8. 



Eqn. ^, (A7Vb(Ty))|G/ 



{ANbiTf))^ ^ mm) 



1 



2A77 



vUT)fiT)dTj26) 



(a) For cj. — 1/3, the fluctuation in net baryon num- 
ber is (A7Vh(T/))2j. - 38.9. (b) For = 0.18, we get 
{ANb{Tf))l^, - 44. (c) Taking from Eqn. |l|, the 
fluctuation in net baryon number is (A7Vb(r/))^Q 
is 158. The evolution of fluctuation for this case is 
shown in Fig. ^. 



V. FLUCTUATIONS AT THE FREEZE-OUT 

The net baryon fluctuation at the freeze-out {Tf 
= 120 MeV) is a combination of the dissipation and 
the generation effect as presented in the previous 
section. The resultant fluctuation is the sum of the 
variances {ANb{Tf))'^ obtained for each of the two 
processes. 

We present results in terms of net baryon fluctu- 
ation per unit baryon, (ANbiTf))"^ / Nb^y. For Pois- 
sonian noise this value should be close to unity. De- 
viation from this numerical values will indicate the 
presence of dynamical fluctuation. 
1. Quark Gluon Plasma: The final fluctuation 
in the net baryon number per unit baryon in the 
QGP scenario at the freeze-out point for the three 
EOS mentioned above are given by, (a) for = 
1/3, {ANb{Tf ))l^^/Nb,y ~ 1.0, (b) for = 0.18, 



{ANb{Tf)Y^^^lNb,y ^ 
Eqn. ^ {ANb{Tf))%^^ 



1.0 and (c) taking from 

/Nb^y ^ 2.0. 



2. Hadron Gas : The net baryon number fluc- 
tuation per unit baryon in the hadronic gas at 
freeze-out for different values of of are, (a) for 



3. Hadron gas with mass variation in medium: 
The net baryon number fluctuation per unit 
baryon for the hadronic initial state with mass 
variation at freeze-out for the three EOS are, (a) 
for c2 = 1/3, {ANbiTf))l^^^,/Nb,y ^ 1.3, (b) for 
c, - 0.18, {ANb{Tf))l^^^,/Nb,y ^ 1.3 and (c) 
taking c2 from Eqn. |l8|, {ANbiTf))^^^^^. / Nb,y - 2.5. 



The net baryon number fluctuations per unit 
baryon is summarized in Table ^ For all the sce- 
narios the values of the quantity [ANb{T})Y / Nb^y 
is larger than Poissonian noise for EOS taken from 
lattice QCD. This indicates presence of dynami- 
cal fluctuations. However, it should be mentioned 
here that the errors in lattice QCD calculations is 
large for temperature below the critical tempera- 
ture. Furthermore (AA^f,(T/))|Q/(AiVf,(T/))|(jp - 

1.4 and (AiVb(r/))2j^,„./(A7Vb(r/))|GP ^ 1-25. 
This means, although scenario (1) may be distin- 
guishable from the scenarios (2) and (3) it is very 
difficult to differentiate (2) and (3) from fluctuation 
measurements. However, for EOS with = 1/3 and 
0.18 they are close to Poisson value of 1. This indi- 
cates that for these EOS; it is difficult to distinguish 
among the three scenarios. 



VI. RESULTS AT RHIC ENERGIES 

The net baryon number fluctuation has been eval- 
uated for RHIC energies of ^/i = 200AGeV of 
Au-|-Au collisions. We have taken the total (charged 
and neutral) dN/dy ^ 1100 The initial tem- 

peratures obtained from the above multiplicity are 
quite large, 370 MeV for hadronic gas and 290 MeV 
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for the case of hadronic gas with mass variation. As 
these temperatures are well above the critical tem- 
perature predicted by lattice QCD, we have consid- 
ered only the QGP scenario for RHIC energies. We 
have taken the initial time to be 0.6 fm/c and 
specific entropy to be 150 for Au + Au collisions. For 
QGP we take gejf = 47.5 which gives Ti ^ 251 
MeV and the constraint on the specific entropy gives 
an initial chemical potential of 73 MeV. With these 
initial conditions for QGP scenario at RHIC, the ini- 
tial fluctuations turn out to be : (A7Vb(Ti))QQp — 42. 
The entropy density in the QGP phase is calculated 
from Eqn. with gefj = 47.5 to be ~ 3900. So 
the fluctuation in net baryon number per entropy is 
{Am{n))lcp/S = 0.011 and {ANt{n))^^^p/m,y = 
1.6. The value of dNh/dy ^ 26 for RHIC energies a 
factor of about 2.4 lower than that for SPS energies. 

As before, here the values of Tc and Tf are 
170 and 120 MeV respectively, (a) For = 1/3, 
the fluctuation in net baryon number due to dis- 
sipation at freeze-out, (AAffc(T/))QQp turns out to 
be 2.7. Whereas the generation mechanism gives 
(AiVb(T/))QQp ~ 35.5. So the resultant fluctu- 
ation in net baryon number per unit baryon is 
(A7Vb(T/))|Gp/iVb,, - 1.46. (b) For cl = 0.18, 
the fluctuation in net baryon number due to dissi- 
pation, (AA^;,(t/))qqp turns out to be 2.3. Whereas 
the generation mechanism gives (A7Vf,(r/))QQp ~ 
48.7. The evolution of the dissipation and the 
generation of fluctuation for the QGP scenario at 
RHIC energies is shown in the Figs. |^ and ||. 
So the resultant fluctuation in net baryon number 
per unit baryon is {ANi,{Tf ))QQp/ Nh_y ^ 1.96. 
(c) If we take as given by Eqn. |l^, the fluc- 
tuation in net baryon number due to dissipation, 
(A7Vb(T/))QQp is 0.25, the generation mechanism 
gives {ANb{Tf))QQp ~ 82.25. So the resultant fluc- 
tuation in net baryon number per unit baryon is 
{Am{rf))l^p/m,y ~ 3.17. 

Note that though the absolute value of the fluctua- 
tions in case of RHIC is smaller than SPS the fluctu- 
ation per baryon at RHIC is larger at the freeze-out 
point. 



VII. SUMMARY 

We have discussed the evolution of the fluctua- 
tion in net baryon number from the initial state to 
the final freeze-out state for three different scenarios, 
viz., (1) formation of QGP (2) hadronic gas and (3) 
hadronic gas with modified mass in the medium. In 
case of QGP formation we have assumed a first order 
phase transition. We find that the fiuctuations at 
the initial stage agree with previously obtained val- 
ues 10 where there are clear distinctions among the 
three cases. The fiuctuations at the freeze-out point 
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FIG. 5: The dissipation of net baryon number fluctua- 
tion for the QGP scenarios as a function of temperature 
calculated for Au + Au collisions at RHIC energies. The 
dashed lines corresponds to results obtained for = 1/3, 
while the solid lines corresponds to results obtained us- 
ing the value of given in Eqn. PLSl 
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FIG. 6: The generation of net baryon number fluctua- 
tion for the QGP scenarios as a function of temperature 
calculated at RHIC energies of 200 AGev Au + Au Col- 
lisions. Notations are same as Fig. ^ 



depend crucially on the equation of state (value of 
Cs). Fluctuations with ideal EOS is seen to dissipate 
at a slower rate compared to EOS from lattice cal- 
culation. At SPS energies the values of the variance 
depend crucially on the EOS. For EOS from lattice 
QCD parametrization the fluctuation at the freeze- 
out is larger than the Poissonian noise. However, for 
ideal gas EOS and EOS with = 0.18, we do not 
observe fluctuations of dynamical origin. At RHIC 
energies the value of the fluctuations are larger than 
the Poissonian noise for all the three EOS under 
consideration here, indicating the fluctuations of dy- 
namical origin. The effects of the finite acceptance 
on the fluctuation enters our calculations through 
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A77 and according to Eqn. and ^ it is same for all 
the three scenarios (1), (2) and (3) discussed above. 
For general discussions on the effects of acceptance 
on the fluctuations we refer to Ref. The de- 

pendence of the fluctuation on the centrality of the 
collisions (impact parameter) get canceled to a large 
extent in the ratio, {ANi,{Tf))QQp/Ni,^y. It is shown 
in 1^ that it is possible to control the impact pa- 
rameter dependence of the fluctuation by measuring 
Et and analyzing data in narrow bins of Et- A 
full (3+1) dimensional expansion will lead to faster 
cooling, and hence it is interesting to see the surviv- 
ability of fluctuations in such a scenario Eq] . 
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